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ABSTRACT 

We present far-infrared spectra and maps of tlie DR21 molecular cloud core between 196 and 671 /im, using the HerschelSVYREi 
spectrometer. Nineteen molecular lines originating from CO, ^'^CO, HCO^ and H2O, plus lines of [N 11] and [CI] were recorded, 
including several transitions not previously detected. The CO lines are excited in warm gas with Tkin ~ 125 K and nnj ~ 7 x 10* 
cm~"^, CO column density N(CO) ~ 3.5 x 10^* cm~^ and a filling factor of ~ 12%, and appear to trace gas associated with an 
outflow. The rotational temperature analysis incorporating observations from ground-based telescopes reveals an additional lower 
excitation CO compoment which has a temperature ~ 78 K and N(CO) ~ 4.5x10^^ cm~^. 
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1. Introduction 

We report observations of the far-IR spectrum of the DR21 
molecular cloud core obtained with the Herschel satellite 
between 196 and 671 /im . The DR21 Hll-region/molecular 
cloud is part of the Cygnus X complex of molecular clouds 
located at a distance of 1.7 kpc (Schneider et al. 2006). 
This region has been subject to numerous studies at dif- 
ferent wavelengths (Richardson et al. 1988, Wilson & 
Mauersberger 1990, Liechti & Walmsley 1997, Schneider 
et al. 2006, 2010, Jakob et al. 2007). The main DR21 cloud 
core has a mass of ~ 20,000 A/q (Richardson et al. 1989), 
and contains one of the most energetic star formation out- 
flows detected, with an outflow mass of ^ 3000 Mq (Garden 
et al. 1991, Cruz-Gonzdlez et al 2007). 



■ ^ 2. SPIRE observations 



X 



2.1. Spectra 

We present Science Demonstration Phase (SDP) obser- 
vations obtained with ESA's Herschel Space Observatory 
(Pilbratt et al. 2010), using the Spectral and Photometric 
Imaging Receiver (SPIRE - GrifRn et al. 2010). The calibra- 
tion and characteristics of SPIRE have been described by 
Swinyard et al. (2010). SPIRE was operated as an imaging 
Fourier Transform Spectrometer (FTS) in the high resolu- 
tion mode (A/AA = 1000 (= 300 km s^^ at 250 ^m) sam- 
pling across an approximately circular field of view with an 
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unvignetted diameter of 2.6'. This means that the line pro- 
files are unresolved. The sky footprint is formed by two de- 
tector arrays: the 19 pixel SLW array (671-303 /Ltm ) and the 
37 pixel SSW array (313-194 /xm ), with beam widths vary- 
ing from 17" at 194 /im to 42" at 671 /im, with uncertainties 
of ± 7-10% (Griffin et al. 2010). The integration time was 
1065 seconds, summed from two seperate observations. The 
current best estimates of the absolute uncertainties for the 
FTS detectors are 10-20% for the SSW detectors, and ~ 
30% for the SLW detectors (Swinyard et al. 2010). 

The unapodised FTS spectra provide the highest spec- 
tral resolution, with a classical instrumental sine function 
line shape. A spectral line fitting routine was developed 
for extracting line parameters (Jones et al. 2009). This fits 
a continuum (either a low order polynomial or a black- 
body variant) using the Levenberg-Marquardt least squares 
method. The fitting procedure weights the spectral inten- 
sity at a given frequency of an averaged spectrum by the 
statistical uncertainty at that frequency, returning line cen- 
ters, intensities, line widths and their associated fit errors. 

2.2. Maps 

The SPIRE observations sparsely sample the field of view, 
although there are calibration uncertainties for the outer 
ring of detectors at the edges of both arrays that are not 
yet fully characterised. To provide a first look at the rel- 
ative distributions in the various species, we have interpo- 
lated the fiuxes of individual lines, although this spatial 
information is not used in subsequent line modeling. 

The maps of selected species are shown in Fig. [2j The 
CO lines in both detector arrays show a prominent cen- 
tral peak, with extensions to the east and west along the 
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Fig. 1. [Upper picture] Spectrum of DR21 with the SLW 
detector array, with beam sizes varying from 42.0" at 671.1 
/im to 37.3" at 316.5 /im, and [Lower picture] with the SSW 
detector array, with beam sizes varying from 18.7" at 281.7 
fim to 16.8" at 194.2 ^m. 

Table 1. Fluxes measured in the central pixel 



Species Transition 



Wave 



Int Intensity 
W m-2 sr-i 



Intens Error 
W m-2 sr"^ 



CO 


J = 4-3 


650.1 


2.85 


-8) 


6.93 


-10) 


CI 


'Pi - 'Po 


609.0 


4.86 


-9) 


9.96 


-10) 


HCO+ 


J = 6-5 


560.5 


3.99 


-9) 


4.29 


-10) 


"CO 


J = 5-4 


544.1 


1.66 


-8) 


5.04 


-10) 


CO 


J = 5-4 


520.3 


6.81 


-8) 


3.39 


-10) 


HCO+ 


J = 7-6 


480.3 


1.02 


-8) 


1.35 


-9) 


"CO 


J = 6-5 


453.5 


2.44 


-8) 


3.21 


-9) 


CO 


J = 6-5 


433.5 


1.15 


-7) 


1.47 


-8) 


HCO+ 


J = 8-7 


420.3 


1.32 


-8) 


2.10 


-9) 


H2O 


2ii-2o2 


398.6 


2.33 


-8) 


3.03 


-9) 


"CO 


J = 7-6 


388.7 


3.66 


-8) 


5.88 


-9) 


CO 


J = 7-6 


371.6 


2.14 


-7) 


1.29 


-9) 


CI 


'P2 - 'Pi 


370.5 


3.03 


-8) 


1.26 


-9) 


"CO 


J = 8-7 


340.1 


6.79 


-8) 


1.80 


-8) 


CO 


J = 8-7 


325.2 


3.15 


-7) 


4.56 


-8) 


CO 


J = 9-8 


289.1 


4.89 


-7) 


4.23 


-9) 


CO 


J = 10 - 9 


260.2 


5.94 


-7) 


1.01 


-8) 


CO 


J = 11-10 


236.6 


7.26 


-7) 


5.46 


-9) 


CO 


J = 12-11 


216.9 


7.44 


-7) 


6.72 


-9) 


Nil 


'Pi - 'Po 


205.2 


1.45 


-7) 


4.71 


-8) 


CO 


J = 13 - 12 


200.3 


6.90 


-7) 


3.96 


-8) 



Fig. 2. The SPIRE maps of selected species toward DR21. 
Since the data points do not fully sample the mapped area, 
this map was constructed by interpolating between the data 
points. The (0,0) position is located at RA (2000) = 20'* 
39™ 01M8, Dec (2000) = +42° 19' 43".66 



well known outflow. This has been assumed to be associ- 
ated with outflowing gas with (Tex ~ 2000 K and H(H2) ~ 
1x10^^ cm~^ from Garden et al. 1991). However, as will be 
seen in the high resolution JCMT observations (Fig. [3]), the 
emission traced in the SPIRE maps is also clearly visible 
in the relatively low excitation CO J = 3-2 data, suggest- 
ing that there may be a mixture of low and high excitation 
gas present. This is confirmed in Fig. [31 where similar ex- 
tensions of the ambient gas are present in the JCMT CO 
J = 3 — 2 map, and that of Schneider et al. (2010). This 
is not unexpected, as this outflow appears to have a very 
large mass of several thousand Mq, and presumably the 
high velocity gas phase overlaps (or may co-exist with) am- 
bient material. The SPIRE maps also show that the ^Pi- 
•^Po atomic carbon line has a similar spatial distribution to 
that of CO. By contrast, the II2O and [N 11] lines appear 
to be more compact and centred close to the DR21 cloud 
core, although the [N 11] distribution is elongated to the 
east - observations with higher signal to noise and better 
sampling are needed for more detailed comparison. 

2.3. JCMT CO J ==3-2 observations 

CO J = 3-2 JCMT archival data (programme M07AU01) 
with a 15" beam and spectral of 0.05 km s~^ are shown in 
Fig. [3l from a 4.5 hour integration using the HARP array 
receiver. The area covered by the SPIRE footprint (Fig. 
[5]) is shown as a white square. The JCMT observations 
clearly trace the outflow which runs from the NE-SW from 
DR21 from the centre of the white box. The JCMT map 
also reveals a prominent north-south ridge that includes CO 
peaks associated with the well-studied sources DR21(0H) 
and DR21-FIR1. Around the DR21 core, a bipolar structure 
close to the systemic velocity is coincident with distribution 
of high velocity gas and shocked H2 (Garden et al. 1991). 

3. Modeling the CO lines 

The most extensive modeling of the CO toward DR21 is 
by Richardson et al. (1986, 1988), Wilson et al. (1990), 
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Fig. 3. CO / = 3-2 integrated emission map of DR21 be- 
tween velocities of -40-±40 km s^^ . The lowest contour is 
drawn at 1 K km s~^ (5cr as observed in a 61 kHz (^-^0.05 
km s~^ ) spectral channel), and the contour intervals are in- 
cremented in steps of 50 K km s~^. The peak antenna tem- 
perature and velocity integrated temperatures at the DR21 
peak are 39.2 K and 726 K km s~^ respectively. The temper- 
ature scale is calibrated in antenna temperature ra*units, 
which is related to the main beam brightness temperature 
Tnibby Tmb = ^a*/ 77nib, whcrc the main beam efficiency, 
'7mb= 0.63. Further details of HARP calibration are given 
in Buckle et al. (2010). 



Schneider et al. (2006) and Jakob et al. (2007). Richardson 
et al. (1986, 1988) presented a multiphase model with gas 
densities spanning the range 10"^-10^ cm"-^, and gas tem- 
peratures in the low temperature component ^ 30 K. Jakob 
et al. (2007) confirmed this using KOSMA and ISO observa- 
tions, finding and additional warm phase component with 
Tkin - 80-150 K and clump density tih^'- lO-^-lO** cm^^. 

Wc initially constructed a rotational temperature dia- 
gram for the SPIRE CO and ^"^CO lines. These were aug- 
mented with the JCMT CO hue from Sect. [231 plus IRAM 
CO J = 2 — 1 observations (Schneider et al. 2010), with 
suitable beam size corrections. The rotational temperature 
diagram is shown in Fig. U) 

Both species show evidence for two gas components, 
a lower temperature phase with a rotational tempera- 
ture Trot = 78 K and total CO column density N(CO) -- 
4.5 X 10^^ cm~^, in addition to a higher temperature compo- 
nent with Trot = 185 K and N(CO) - 9.7x10^^ cm'^. The 
^■^CO lines are more limited and noisy, with the SPIRE 
lines indicating an intermediate temperature phase having 
T,ot= 109 K and N(i3C0) ~ 8.8x10^^ cm'^. The data 
for "'^^CO also show evidence for a low temperature compo- 
nent, although this relies on comparison with low frequency 
ground based data (JCMT, IRAM) obtained with different 



1 T„, = 78 K N„| = 4,5 x 10" crr 
\ CO 



^f--,. 



t--. 



T,„, = 185 K N„, = 9.7 x 10" cm=: 



^'CO J- 



T,„- = 109 K N^„, = S-8 X IC err 



'00 200 300 400 500 

Upper state energy level E,^, (K) 

Fig. 4. CO and ^"^CO line rotational temperature diagram 

beam sizes. Such a result is expected, since the observations 
probe deeper into the PDR of each clump in ^■^CO than in 
CO. 

There are several problems with the rotational temper- 
ature approach, including wavelength dependent beam size 
corrections, opacity and calibration errors. These uncer- 
tainties can however be mitigated by i) taking ratios of 
the various CO line intensities on a single detector and us- 
ing these to constrain the excitation conditions though our 
LVG modeling, and ii) using observations from the central 
pixel where the SSW and SLW beams are coincident and 
the calibration is well determined. This approach particu- 
larly mitigates against the beam size and calibration errors, 
since only fiux ratios are being used to estimate the excita- 
tion conditions. 

The model fit was made to the CO and ^^CO lines 
using the off-line version of the RADEX LVG code (Van 
der Tak et al. 2007). The line ratios observed on the same 
detectors (hence the beam sizes are similar) were used to 
restrict the likely excitation conditions. It proved difficult 
to find an unique single temperature model that simulta- 
neously predicted the relative intensities of both isotopo- 
logues. However, the SPIRE data can be approximately re- 
produced by a single phase moderate temperature gas with 
Tkin- 125 K, volume density ~ 7 x lO^* cm^^^ with N(CO) 

- 3.5 x IQi^ cm-2, filling factor ~ 12%, and a [C]/[i3C] ra- 
tio of 65. This model does however slightly overpredict the 
low- J ^'^CO 4-3 - 6-5 line intensities, compared to the J = 
7-6 line. Changing the temperature and density from these 
conditions considerably worsened the high- J CO line fits, 
although a more complex multiphase model, with appropri- 
ate (and uncertain) beam size corrections would improve 
the fit of the low J-lines. We have not attempted to fit to 
a PDR-model, as the data and calibration quality need to 
be improved if tests between models are to be made, and 
that this is beyond this first look paper. 

4. Modeling the H2O and [N 11] lines 

An objective of this study was to detect the [N 11] 205 /im 
line, and to compare it with the [C 11] 157 line which has vicrit 

— 46 cm^'^, Te — 8000 K. This has a nearly identical criti- 
cal density for excitation in ionised regions. Their line ratio 
is directly related to the N+/C~'' abundance ratio, and this 
ratio traces the fraction of the observed [C 11] emission that 
arises from ionized regions (Oberst et al. 2006). Taking the 
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Fig. 5. LVG analysis of the CO and ^'^CO lines toward 
DR21. The SPIRE data points, (J = 4-3 to 13-12) are 
shown as filled circles, and the ISO data points (Jakob et 
al. 2007) as squares, along with 30% absolute error bars. 
The CO J == 3-2 fine was obtained from the JCMT data, 
convolved to a resolution of 17" to match the beamsizc of 
the nearby SPIRE lines. The filled circles are the SPIRE 
data, open sqares are ISO data from Jakob et al. (2007), 
open triangle is determined from the JCMT data in this 
paper, and the crosses are the SPIRE ^"^CO data. The top 
solid curve shows the fit to the predicted SPIRE CO line 
intensities, and the lower dashed one is for ^■^CO with an 
assumed CO/^^CO abundance of 67, an estimated filling 
factor of 12%. As might be expected, the fit for a single 
tempeature lies intermediate between the two Trot values 
inferred from Fig. U) 



SPIRE upper limit of 7.5 x 10"^ W m'^ sr"! with Jakob 
et al. (2007), the ratio of the 122/205 nm lines is > 1.9, 
which is only adequate to constrain the ionised gas density 
to be > ^ 30 cm~^. The [C ii] /[N ii]205 ratio using the 
Jakob et al. (2007) tabulation is > 5.6. Given current un- 
certainties and lack of an [N Ii]i22 flux, it is necessary to 
await improved data. We note that the [N ii] extension to 
the east (see Fig. H]) coincides with a hole in the excited II2 
emission image (Cruz-Gonzalez et al. 2010), which may in- 
dicate there is a cavity of ionised gas. However clarification 
will require future observations with better sampling. 

In Fig. [6] we show a section of the spectrum with the 
398.6 /xm para-H20 line, and the HCO+ J = 6 - 5, 7 - 
6 and 8 — 7 lines. Putting the SPIRE sensitivity into per- 
spective, Jakob et al. (2007) report that the integrated CI 
'^Pi-'^Po intensity measured from the KOSMA telescope 
with an 80" beam is 46.6 K km s~^, and main beam bright- 
ness temperature '^ 25 K. By comparison the same line ob- 
served with SPIRE has a peak S/N ratio of ^14 as seen in 
a single SPIRE channel. We also used RADEX to compute 
an LVG solution for the 2ii-2o2 para-H20 line at 398.5 /xm. 
Assuming similar excitation to that from the CO solution, 
for an abundance X[Il20] = ^ 4 x 10~^ and line width 
of 40 km s^^ (Hjalmarson et al. 2003), we predict that the 
SPIRE fiux should be 2.4 x 10"^ W m"^ sr~\ which agrees 
with the measured value of 2.33± 0.3 x 10"^ W m"^ gj.-! 

5. Conclusions 

We have presented the SPIRE spectrum of a star-forming 
molecular core, DR21, showing for the first time the com- 
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Fig. 6. Continuum subtracted spectrum showing the 398.6 
fira para-Il20, and HCO+ J = 7 — 6 and J = 8 — 7 lines. 
Some instrumental features remain - notably a broad bump 
close to the CH+ line at 347.9 /im. Although it is possible 
to fiatten this spectrum further, we chose only to fit a third 
order polynomial to the entire 671 -319 /im range of the 
SLW detector. 



plcte CO and ^^CO band head from J = 4-3 to 13-12, along 
with their maps at far-infrared wavelengths. A rotational 
temperature analysis shows two gas phases with Tj-ot ~ 80 
K and CO column density ~ 4.5 x 10^® cm~^, and Trot = 
185 K and N(CO) --10^^ cm~2 respectively, although the 
-'^'^CO Trot is intermediate between these two. Simple LVG 
modeling shows the presence of warm (125 K) and dense 



(nH2~ 7 X 10 cm""*) gas, which is traced by in the SW 
extension in the CO and CI maps. The observed flux from 
the 398.6 /im II2O line is consistent with these values. 
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